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Abstract 

We show how to define and go from the spin-s spherical harmonics to 
the tensorial spin-s harmonics. These quantities, which are functions on the 
sphere taking values as Euclidean tensors, turn out to be extremely useful 
for many calculations in General Relativity. In the calculations, products 
of these functions, with their needed decompositions which are given here, 
often arise naturally. 



1 



I. Introduction 



The use of the ordinary spherical harmonics, Yi m (8,(p) (and their as- 
sociated vector and tensor harmonics) with all their properties, their eigen- 
value/eigenvector behavior, their orthonormality properties, their use in solv- 
ing various problems and equations in mathematical physics, their interpre- 
tations in terms of multipole expansion have been ubiquitous in theoretical 
physics. Several years ago a generalization of the ordinary harmonics was 
developed^ and referred to as the spin-s spherical harmonics and denoted 
by 

(s)Yl m {0, V 3 ) 

with s = the ordinary spherical harmonics. They have proved to be very 
useful in problems involving spin-s fields and have become an almost essential 
tool in problems involving gravitational physics. Though the ( s )Yi m are closely 
related to the Wigner D-matrices, D™" 1 {and to generalized versions of vector 
and tensor spherical harmonics} and can be derived from them, nevertheless 
it is the particular form and specific properties of the ( s )Yi m that have been 
of great use. 

Recently, it turns out, that there are quantities that have been appearing 
in many calculations that are closely related to the ( s )Yj m , which we refer 

(s) 

to as the tensor spin-s harmonics and are denoted by YLi k . This note is 
devoted to a discussion of these new harmonics. 

(s) 

There is a one-to-one correspondence between Y^J k and the ( s )Yj m . 

The indices i k indicate symmetric and trace- free 3-dimensional Euclidean 

tensors; the number of tensor indices equals I. The number of independent 
components, in both Yk] k and ( s )Yj m , is the same; N = 21 + 1 and, in fact, 

(s) 

the quantities Y^J k are just linear combinations of the ( s )Yj m and could 
in principle be written as 

k — ^mJ^^i k (s) * tra- 
it is however difficult and a bit unwieldy to find the ' s and a more direct 
approach, namely to define them independently, is considerably easier. 

Is) 

The value in using the YhJ k instead of the spin-s spherical harmonics 
is two-fold, (i) They simply appear as they are defined in many of the calcu- 
lations and (ii) in non-linear theories very often within calculations one finds 
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that products of these tensor harmonics automatically are there. They, in 
turn, must be decomposed via a Clebsch-Gordon expansion. These expan- 
sions are most easily done with the Ykj k . 

In section II we will describe how the vif] k arise in GR and precisely 
how they are defined. In addition for clarity we give many examples. In 
section III, we discuss how products [Clebsch-Gordon expansions] are found 
and give an example from GR in section IV. The main bulk of this note 
appears in appendices. Appendix A contains some useful miscellaneous re- 
lations while Appendix B contains a table of Clebsch-Gordon expansions of 
some of the most useful products. 



II. The Tensor Spin-s Harmonics 



A. Notation & Basics 

We begin with either Minkowski space or at an arbitrary point in a 
Lorentzian space-time and take the metric as 

r] ab = diag[l, -1,-1,-1]. (1) 

The standard null tetrad, parametrized by the complex stereographic 
angles (£, (), with ( = e 1 ^ cot |, can be given by 



l a = ^ (i + CC, C + C, -»(C - 0, -i + CC) , (2) 

^^(0,1-^(1^), (3) 

^^(0,1-^(1 + ^), (4) 

7I(ITc7) ^ + ~ (c + ^' ^ - ^' 1 - ' (5) 



m a 



m a 



n a 



As (£, C) sweeps out the sphere, the null vector l a sweeps out the light-cone 
and 'drags along with it' the remainder of the tetrad. From the tetrad we 
obtain the time-like and space-like vectors, 
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t a = l a + n a = ^(1,0,0,0) (6) 
c a = l a -n a =-^-(0,C + C,-i(C-C),-i + CC) (7) 

with norms 2 and -2 respectively. 

The main tool or ingredient in our work will be the 3-dimensional {Eu- 

clideanjvectors, = 1, 2, 3], obtained by projections that are normal to 

r, i.e., 



' =-(c + C,-*(C-C),-i + CC), (8) 



v^i + CC) 

mi = ^^(l-f^l + cVc), (9) 

m = -^L_(i-e,i(i + a2c), (10) 

m = -k= TfTf^j ((C + C),-»(C-C),-i + CC), (11) 

c i = ii-n< = ^^(C+?,-<(C-?),-l + Cf)- (12) 
In terms of # and 0, 

c, = — v^2(cos 0sin 6, sin 0sin 6*, cos 6). (13) 

and hence q is just — \/2 times the unit Euclidean radial vector. Note that 
we have used the Minkowski metric, (JH to raise and lower even the Euclidean 
indices in Eqs.(??) - (|T2l which gives rise to the minus sign. 

Some of the important algebraic properties of these vectors obtained by 
direct calculation from their definitions, are 

1 _ _ 

dij = —CiCj + mirrij + rriimj , (14) 

c fc = -V^iekijiUiMj, (15) 
m k = -^ekijrriiCj, (16) 
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% 

nijm k - rrijink = -j=e jki Ci, (17) 



m k Cj - rrijCk 



iV2e kji mi. (18) 



Noting that Ck,m k and ra k have respectively spin-weights (0,1, -1) and an 
/-value, I = 1, we have, either from the definitions of the differential operator 
edth acting on a spin-wt. s function rjr s \, i.e., 

h V(s) = Pt s d c (P^ is) ), (19) 

h V(s) = Pt +s %(P s V (s) ), (20) 

= 1 + CC, (21) 
or from the general eigenvalue relation^ 

55 (s) lk = -{I - s)(l + s + l) (s) Y lm , (22) 
95 (s) ^ m = -(l + s)(l-s + l) (s) Y lm , (23) 

the differential relations 



5q 


= 2mj, 


(24) 




= 2m j, 


(25) 






(26) 


5mj 


Q) 


(27) 


55q 


= -2c, , 


(28) 


55mj 


= -2m;, 


(29) 


55mj 


= -2m;. 


(30) 



B. Definition of Y, 



The essential idea to define the tensor harmonics YnJ k is simply to use, 
in an appropriate way, tensor products of the three basic Euclidean vectors, 
(q ,mj,m"j). Since each of them have the value I — 1, in any product of n 
terms the /-value of the product is / = n. The spin wt. of the product is given 
by the algebraic sum of the spin wts. of the constituent vectors. In addition 
we require the products to be symmetric trace-free Euclidean tensors. 
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The easiest way to do this is to first define the spin-wt. s = I tensor 
monic k , {with I indices, i....k or / factors m,k] and Yk^_ k by 



Y (l)i k = m i m 3 m k (31) 

Y (i)iL...k = rn i m j ...m k . (32) 

First note that both are obviously symmetric and trace-free [since rriimjS 13 
0] and that any derivatives will also be symmetric and trace-free. Further- 
more we recall that the edth operators 

5 and 5 



are stepping operators, i.e., they add or subtract one to the value of s. 

A* 
(0 



For positive values of s, [s = 0,1, ....,/], the Ykj k are defined by ap- 



plying the operator 

5, 

/ - s times to k , i.e., 



k = Z l - S {Y$ t J (33) 



and for negative values of s, i.e., for [0, —1, —I] 

(V)i k ~ ° I 1 (I) 



Y&iz* = m^;; :':....,}• (34) 



Note that 

It is easy to show that the tensor harmonics satisfy the eigenvalue equations 



v i-\s\) _ -y (s) 

k — 1 k- 



ZZY$ k = -(l-s)(l + s + l)Y$ k , (35) 

53y ( Sj; k = -(l + s)(l-8 + l)Y$ k . (36) 

C. Examples 

In order to clarify these definitions and for later use, using Eqs. (|2^|) - (|27j) 
and {TJJ, we give several examples: 
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r — i. 



Y 1 
Y° 



mi, 
mi. 



Y£j = rmmj, 

Y iij = Sy 2 2 y = -{ciTUj + rrncj), 

Y^ij = dY 2i j = 3ciCj — 25ij, 

Y 2 ~j = -(ciMj+miCj), 

Y 2 ~f = miWy 



miirbjm k 

—(cirrijmk + rriiCj-mk + rriimjCk) 

5 

-[Sijirik + Skjirii + Sikirij] + - {ciCjUik + c k Cj 
6(5ijC k + S ik Cj + S kj Ci) - 15ciCjC k 



Y tijki = rrumjmkmi (47) 

Y 4ijki = -Cimjm k mi - miCjm k mi - mimjC k mi - mimjm k Ci 

Y^iju — 2[ciCimjm k — miminijm k — m k mimjmi — rrTjmim k mi — mimjm k mi 

+c j m k {cim i + Cirrii) + c k {cim(mj + (cjuii + Cim^mi)} 
Ylijki — Q[mic k mimj + m k cimimj + mjCimim k + miCimjm k + mjC k mimi 

+miC k mjmi + Ci(—c k Cirrij + mirrijm k + (rn k rrij + mjm k )mi) 

+Cj(-CiCim k + Mimim k + m k mimi + mim k mi - c k (cirrii + Qm;))] 
Y 4ijki = l05ciCjC k ci + I2[5 ik 5ji + 5u5 jk + 8^8^] 

-30[5ijC k ci + S ik CjCi + Suc k Cj + S kj CiCi + 5ijC k Ci + 5 k iCiCj] 



III. Products of the 



Very often in doing detailed calculations in general relativity due to the 
non-linearity one has to deal with products of different Ykj k . Usually they 
involve only small values of both s and /, most often from to 3 or 4. In 
principle, by using the products of the Wigner D-functions, i.e., Clebsch- 
Gordon expansions, one could work out the Ykj k products. In practice 

(s) I 

since the conversion of the Y„J k to the D mm , is quite complicated - with 
various conventions - and we are most often interested in the low values of s 
and Z, it is easier to do the calculation of the expansion directly. 
First, instead of using (c* , m^m^), we use 

(Y 1 ° l ,Y 1 1 i Y u 1 ) = (-c l ,m l ,m l ). 

The simplest and most important product relations are found most easily 
by direct calculations using the definitions, (|9|. (|T0l) and (JT2l) . of (q ,mj,mj) 
yielding 



^li^lj 1 ~ Y lj Y U l ~ ~~^ e ijk Y i k , (48) 

Y^-Y^Y* = iV2e jki Yl (49) 
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To find the products Y^ k Y^, YJ^Yy and Y^Yy one writes them out as a sum 
of tensor terms from I = to I = 2 and applies the operator 



several times using the eigenvalue equations (j35|) and |36|) to determine the 
coefficients. 

For example Yy Y^ could be written out as 

= AY ° + + &%% (50) 

with A, B 1 and C 4J to be determined, and then from the eigenfunction rela- 
tions, we have 

55^°] = -2^Y°-6^ fe Y 2 % ( 51 ) 
mm\Y^ k \ = AB% j + 36C j %%. (52) 

Since the left sides of (jSDJ , (JHIJ and l|H2J) are known we can evaluate the A, B l 
and C lJ yielding 



Y^ = -5i 3 + -Y« l3 . (53) 



2 r 1 

3 lJ 3 

In a similar manner, with the help of ijlHJ) and ijlHJ) . we obtain 



^li^ij 1 — ~^~ e ijkYi k Y^^2ij- (55) 



In principle all products expansions can be found in this manner. As the 
details can become quite tedious we simply list the most important ones in 
the Appendix B. 



IV. Application: The Robinson- Trautman Equation 

Though there are many applications of these results that will be given 
elsewhere, here we will show how the tensor harmonics can be used to ap- 
proximate solutions to the Robinson- Trautman equation,^ the final equation 
that determines the type II, non-twisting metrics. Our purpose, in this ex- 
ample, is simply to illustrate how the tensor harmonics and their products 
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enter into GR calculations and so we will not be concerned with the certain 
details. 

I. Robinson and A. Trautman in their investigation of algebraically spe- 
cial vacuum metrics with non-twisting principle null vector found that the 
algebraically special type II metrics could be reduced to the single partial 
differential equation for a 'mass parameter', x(r), as a function of the 'time' 
parameter and a function 

P(r,CX) = v(t,c,C)(i + CC) = V(t,CX)Po, 
a time dependent conformal factor for a two-surface metric. 

Remark 1 The function x( T ) though not the Bondi mass is a close relative. 
The Bondi mass is given, up to a numerical factor, by Mb = x( r )W^( r )i 
where _ 

One can easily show* from Eg. \57a\) that M' B < 0, i.e., the Bondi mass 
loss theorem. 

The Robinson/Trautman equation could be rewritten,^ using the edth 
notation, (fTHl) . as 

X - 3yX = V 3 [SmW + 25"5V] - y 2 (aV)(5"V) (57a) 

with ( ' ) meaning the r-derivative. They showed that by choosing a different 
time parameter t —>■ r* = F(t) and rescaling the V one could then make 

x' = o 

thereby simplifying the equation. We however will use the reparametrizaion 
freedom in a different way: we choose it to make the leading term in V to 
be one and thereby allow x to be time dependent. In particular we assume 
that V takes the form of the tensor harmonic expansion with the I = term 
unity: 

V=l~^Y 1 ° l + ^Y 2 % + (58) 

with rfi symmetric and trace-free. 
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Remark 2 ; In other publications, and for specific reasons, we have chosen 
the T) 's as r -derivatives of other functions, i.e., r/ = In the present work 
this is not necessary. 

Our object is to try to find approximate solutions to (|57a|) by assuming 
that, in an expansion near the Schwarzschild solution, \ is zero order while 
rf and rfi are respectively first and second order. Note that the order of the 
T-derivatives is at this stage not know but be determined by the differential 
equation, (|57aj) . Furthermore we will truncate the harmonic series at two, 
the quadrupole term. The expression for V, (|58|) . will be substituted into 
(I57a,l ) and expanded up to the 1 = 2. This yields three different evolution 
equations, [I = 0,1,2], for the r-derivatives of x, rf and rfK For a reason 
made clear later we will keep terms up to fourth order, even though most 
will eventually be discarded. 

As a preliminary to the substitution of V into (|57ajl we calculate from 
the eigenvalue equations, A35I) and ijHHJ) . the following relations: 

b W = —rf 5oY° + V mY 2 % = rf F° - Q^Y 2 %, (59) 

mW = -2rf Yg + 367/^4 , (60) 

5555V + 2 55F = 2Arf i Y^ ij , (61) 

tfV = 5 2 [1 - ~rf Y-x + rpY*,] = V ^%% = 2A^Y 2 %, (62) 

§V = 24^Y 2 T 2 , (63) 

(5V)(5V) = {2Afrf^ kl Y 2 %Y 2 ll (64) 
Using these relations, Eq. (|57ap can be written as 

V X ' - W'x - V 4 [2A^Y 2 % + ...} + V%2Afrf^ kl Y 2 %Y 2k ^ + ...] = (65a) 
By expanding V 4 as 
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V 4 = [l-i^ + ^l* (66) 
= 1 - 2rfY* + \ifr?Y°Y° + A^Y 2 % + ... 
and using the product 

V&ii = + \y 2 % (67) 

we have 

V 4 = 1 - 2t/* + 77*7/* + rf + H y )YL + - (68) 
Up to fourth-order, Eq. ()65a|) reduces to 

-X + ~(xV - 3^')n° + l-XV ij + 3xV + 24(1 + r/V )7fK(69) 

-W^y^ + (1277V' + seTfV'y^ 
-(24) W^ 2 ^ 2 = 

From the product relations 

y^jk = -pkjY& + l(S tJ Y? k + 5 ik Y^) + l -Y» jk 

and those for Y 2i ^Y 2kl and Y 2 2 i jY 2k f, given in the appendix B, (and using the 
symmetry and trace-free properties of rf- 7 ), we have after a lengthy calculation 
that Eq. ([69|) becomes, to fourth-order and to the I = 2 harmonic, 

-X' + ~{xV - 3xr/*'}y° + + 3 X ^ 7 + 24(1 + r/VV^IjO) 

+48{jjfaV + - f {[^VW - (~r/V + ^V^IX?!) 

_4 8{ ^i^yo } _ (24) 2 {1^^ - ^7/*y fe y 2 y = o. 

From the I = 0, 1, 2 harmonics this is equivalent to 
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-X+^(\vW +3v ij )v ij = 0, (72) 
2(24) 2 

3xrf' ~ xV + -^-rfrf 1 = 0, (73) 



-x'v iJ + 3rf' + 24^' + r? fc rf j (74) 

+SymTrFr(ij) {^j-{^rfr] k + brf k )r] k ^ = 0, 

where SymTrFr(ij) means symmetrize and remove the trace. 

These three equations will be used in two different ways: (i) They will 
first be used to determine the order of the r-derivatives and consequently 
simplify the equations and (ii) then integrate them. 

From the first we see that \ is fourth order 

X=^Y~{\rfrf +3rf)^ = 0. (75) 

Using this in the second equation and keeping only the leading term, we have 
that rf' is third order, 

(76) 

while the third equation for rf*', is second order, 

xr ff = -tofi. (77) 

These equations, by taking them in the reverse order and treating \ — con- 
stant in the 2 nd and 3 rd equations, can be integrated as 

rfj = rfie- Sx ~ lT . (78) 

The second equation can be integrated exactly with the solutions depending 
on an exponential of the exponential e _8x T , i.e., it involves functions of the 
form 

cxp[-8x~V] 
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which goes to one rapidly so here we simply take rf as constant, 

rf = Vo- 

[The same result could be obtained by neglecting the third order term.] The 
first equation then simply integrates to the form 

X = Xo + Ke- m ^\ (79) 

Since here we are not interested in the details of the Robinson-Trautman 
metric, but only to illustrate the use of the tensor harmonics, we will not 
give the detailed expressions for these quantities. We simply mention that 
using Eqs. (|56l) and (|79|) . we can obtain the evolution of the Bondi mass. 
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VI. Appendices 

A. miscellaneous relations 



(1962). 



1: 



(80) 
(81) 
(82) 
(83) 
(84) 
(85) 
(86) 



ii ii 

= ^551^ = 

5K7 1 = 55K7 1 = o 



ay° = -2rm = -2y^ 

_SY° = -2m i = -2Yu 1 
88Y* = -2F° 
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my* = 
mr- 1 = 



—2Y 1 
—lY^ 1 



(87) 
(88) 



Yxi = W 2% = -{c l m J + m l c j ) (89) 

W 2 % = 881* = -4m,m, = -AY 2 % (90) 

Y&S = ^ Y h = 3c ^ - 2 *« (91) 

W 2 % = 6{m i c j + c i m j ) = -6Y^ (92) 

88Y° = -6Y° (93) 



1 = 3: 



Y 3ijk = d Y £ijk = -(cimjiiik + rriiCjmk + m^Ck) (94) 
Y lijk = ^ Y 3%k = 2[ciCjm k + dTUjCk + rriiCjCk 

—rfiimjmk — mim^my, — mimjTn k } (95) 

Y 3ijk = ^ij m k + Skjmi + SikTHj] 

— 10[mimjm,k + mimjmk + mimjfnk] (96) 

5 

Y^ ijk = -[5ijm k + 5 kj mi + 5 ik mj] + -[ciCjm k + c fc c,mj + CiC k mj](97) 

Y 3%k = ^ Y iijk = 6($ijCk + SikCj + tfjyCj) - 15CiCjC k (98) 

Wi ijk = MY£ jk = -Q mi m 3 m k = -QY 3 % k (99) 

<SY 3 ) jk = mY* jk = -10Yi ijk (100) 

813,-* = -12Yi ijk (101) 

881** = -10Yi ijk (102) 

88K^ = -10Yi ijk (103) 

88Y* fc = -12Y 3 ° ijk (104) 
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B. Clebsch- Gordon Expansions 

We present a table of products involving functions with s = (2,1,0, —1, —2) 
and / = (0,1,2). 

1. Products of I = 1 with 1 = 1: 



Y&Zt = ^ l]k Y^ k + l -Y^ (105) 

Y^f = \k 3 -^^-\-Y 2 % (106) 

*M- = l 5 v + l Y l ( 107 ) 
Products of I = 1 with 1 = 2: 



Y&% k = Y 3 % k (108) 

Y£Y 2 % = -^^ + 1(6^ + 6*^) + ^ (109) 

9 3 3 

ylyO _ _ V X A _ V 1 A 

I U I 2jk — r ) I li°kj ^ I lk°ij 

+^ikiY 2 ) l + e tjl Yi kl ) + ^ (110) 
YMjk = -\WuY^ k ] (HI) 

V-iv 1 — V°A 4- V°A _ 2 V°A 

I 2ij I ik — iQ u i k |q ir ifc 10 J ' 

+ e^-] - (112) 

OS* = + + 

i - 4 l 

_ ^7| ( e ikiY 2 ji + £ijiY 2k i) + Y^^jfe (113) 
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Y 2 %Y? k = m^ik] (115) 



3. Products oil = 2 with I = 2 
total s = 4 



total s = 3 



2 1 * 3 3 1 3 

total s = 2 

The numbers in parentheses as superscripts give the eigenvalues of the asso- 
ciated quantities: 



yi y i 3 2(0) 4 2(14) 

I 2kl I 2ij — rj^ijkl "I" rj^ijkl 1 

yo y 2 3 /r 2(0) 4. 1 /r 2(6) 4 3 /r 2{14) 

I 2kl I 2ij — rj-^ijkl 2 V kl rj^ijkl > 

where the K's are eigenfunctions with Z = (2,3,4) 

2/yO y2 , y2 yO \ 
^\ I 2kl I 2ij ' I 2kl I 2ij)' 

I 2kl I 2iji 

I/yO y2 , y2 yO \ 
\ I 2kl I 2ij ^ I 2kl I 2ij)i 



^ijkl 


— I 2kl I 2ij 


^ijki 


_ y-0 V 2 
_ I 2kl 1 2ij 


^2(14) 
^ijkl 


— Y 1 Y 1 4 

— 1 2kl 1 2ij + 





(118) 




(119) 


/ = 2 


(120) 


/ = 3 


(121) 


Z = 4 


(122) 



and explicitly decomposed as 
4$ = -j[(^ + ^^ 



^2(6) _ J_r y2 



^Jfci ^ — Q^4ijW) (123) 
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total s — 1 



yO yl _ 3 xl(0) 1 ,1(4) 1 .1(10) 3 i(i 8) 

I 2ij I 2kl — ^Q J 4ijkl "r -y^Aijkl ~r ^ J Aijkl ~t~ j J 4ijkl K 1 ^) 

y-ly2 1 rl(0) 1 ,1(4) _ 1 .1(10) _ 1 .1(18) n 9 ^ 

I 2ij I 2ki — 20 4ijH 28 3 ' M 20 4ij ' fci 2g 4i i kl y 1 ^ ) 



where the J's are eigenfunctions with / = (1,2,3,4) 



with 



Tl(0) 
J 4ijkl 


r 1 


Jy 4ijkl 


+ 4-Eli jM - 4F 4 1 ijH , 


/ = 


1 


t1(4) 
J 4ijkl 


r 1 

— ^Aijkl 




+ l ^ E lijki + ^IFlijki > 


/ = 


2 


.1(10) 
J 4ijkl 


r 1 

— ^Aijkl 


Jy 4ijkl 


_ ^Xijki + S-Fijfcj , 


/ = 


3 


.1(18) 
J 4ijkl 






_ ZElijki ~ ^lijki i 


Z = 


4 



r 1 

^njki 


= yO yl 

— I 2ij I 2kl 


n 1 

Jy 4ijkl 


= yi yo 

— I 2ij I 2kl 


TP 1 

nj 4ijkl 


= Y^Y 2 

— 1 2ij I 2kl 


TP 1 

r 4ijkl 


— 1 2ij 1 2kl 



(126) 
(127) 
(128) 
(129) 



(130) 
(131) 
(132) 
(133) 



and decomposed as 



ifc e «j/ + djk e lif + $il e kjf + ^7 e fei/]^l/) 

.1(10) _ _ 1 r V l V l V l , yl 1 

J 4ijfcZ — [ e jfce-^3 ej 7 "I" t jke I 3eil t Ue 1 3ejk t jle I 3eik\i 

total s = 



V 2 V- 2 - 1 _l_ 2 P°( 2 ) _L 2 P°( 6 ) 1 P 0(12) 1 p 0(20) 

1 2kl I 2ij — ^ r ijkl "T p^jjfcZ "T ^ijfci "T ^Q P ijkl jQ r ijkl ' l i ' 30 J 
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I 2kl 1 2ij — ^ r ijkl "I" g^r/fci y^ijfci ejijfc* 35 ijM ' ^ ' 



I 2kl I 2ij — g ^ijjfci 7 r jjJW "I" 35 «jH ' l i<5 'J 

with the eigenfunctions F with / = (0, 1, 2, 3, 4) 



: ^ijiM — (^2fe/^2i/ + Y 2ij Y 2kl) + . (^2fe/^2i/ + ^2iy^2u) + 0/) ^2fe/^2^38) 



1 : -^jjw = (^2W^2i/ ~~ Y 2% Y 2kl) + g(^2fe/^2i/ ~~ Y 2ij Y 2k])i (139) 

2 : ^i/fe/ = ( Y 2kl Y 2ij + ^2ij^2feJ 2 ) ~~ g(^2M^2i/ + Y 2ij Y 2kl) ~ ^^2W^2^40) 

3 : ^i/fej — (^2W^2i/ ~~ ^2ij^2fcf ) ~~ ^(^M^W ~~ Y 2ij Y 2kl)i (141) 

4 : -^j/fe/ — O^ikl^ij + ^2ij^2fcf) ~~ 0^2kl^2ij + Y 2ij Y 2kl) + 7^fcZ^2ij)(l42) 



and decomposed into 



Fijkl = i—^(3jl e ike + 5ki£jle + Cjke^U + 5kjeile)Y± e , 

F ijkl = Q^ij Y 2°kl + bklYiij) - g(^j^2ifc + ^fci^lj + $li Y 2°kj + ^fcj^2ij)) 
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